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ABSTRACT 



We analyze the behavior of the parsec-scale jet of the quasar 3C 454.3 during pro- 
nounced flaring activity in 2005-2008. Three major disturbances propagated down the 
jet along different trajectories with Lorentz factors P >10. The disturbances show a 
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clear connection with millimeter-wave outbursts, in 2005 May/June, 2007 July, and 
2007 December. High-amplitude optical events in the i2-band light curve precede peaks 
of the millimeter- wave outbursts by 15-50 days. Each optical outburst is accompanied 
by an increase in X-ray activity. We associate the optical outbursts with propagation 
of the superluminal knots and derive the location of sites of energy dissipation in the 
form of radiation. The most prominent and long-lasting of these, in 2005 May, oc- 
curred closer to the black hole, while the outbursts with a shorter duration in 2005 
Autumn and in 2007 might be connected with the passage of a disturbance through the 
millimeter- wave core of the jet. The optical outbursts, which coincide with the passage 
of superluminal radio knots through the core, are accompanied by systematic rotation of 
the position angle of optical linear polarization. Such rotation appears to be a common 
feature during the early stages of flares in blazars. We find correlations between optical 
variations and those at X-ray and 7-ray energies. We conclude that the emergence of 
a superluminal knot from the core yields a series of optical and high-energy outbursts, 
and that the mm- wave core lies at the end of the jet's acceleration and collimation 
zone. We infer that the X-ray emission is produced via inverse Compton scattering by 
relativistic electrons of photons both from within the jet (synchrotron self-Compton) 
and external to the jet (external Compton, or EC); which one dominates depends on 
the physical parameters of the jet. A broken power-law model of the 7-ray spectrum 
reflects a steepening of the synchrotron emission spectrum from near-IR to soft UV 
wavelengths. We propose that the 7-ray emission is dominated by the EC mechanism, 
with the sheath of the jet supplying seed photons for 7-ray events that occur near the 
mm- wave core. 

Subject headings: galaxies: active — galaxies: quasars: individual (3C 454.3): galaxies: 
jet — polarization: X-rays: galaxies 



Introduction 



During the past four years, the quasar 3C 454.3 (z=0.859) has displayed pronounced vari- 
ability at all wavelengths. In spring 2005 it returned to the night sky with unprecedented bright- 
ness, R ^12.0 mag , a level not seen at optical wavelengths over at least 50 years of observations 



(jVillata et al.ll2006l ). An increase in activity occurred at X-ray and radio wavelengths as well, with 
the 230 GHz radio v ariations having a delay of ~2 months with respect to the optical variability 
2008bl ) . This prominent outburst was followed by a more quiescent period at all 



(Raiteri et al 



wavebands from spring 2006 to spring 2007. During this interval, the optical spectrum possessed 
characteristics typical of radio-quiet active galactic nuclei (AGN), such as a "big blue bump" and 
"little blue bump," attributed to thermal e mission from the a ccretion disk and broad emission 
lines from surrounding clouds, respectively ([Raiteri et al.l 120071). After th e quiescent state, the 
quasar underwent a new stage of high optical activity (iRaiteri et al.ll2008al ) that continued to the 
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end o f 2008. During this time span, very bright 7-ray emission was detected (IVercehone et al 
20081 : iTosti et al.l 120081). with a n excehent correlation between the 7-ray and near- infrared/optical 
variations (jBonning et al.ll2009l ). 



Models proposed to explain the observed variability and spectral energy distribution (SED) of 
3C 454. 3 across the e l ectrom agnetic spectrum involve processes originating in the radio jet of the 
quasar. IVillata et al.l (j2007l ) suggest that the very high optical flux in spring 2005 was connected 



with a disturbance (e.g., a shock) propagating along a curved trajectory in the jet, with optical 
synchrotron photons emitted over a different volume than the longer-wavelength radiation. As 
the emission zone of a given wavelength passes closest to the line of sight, the flux peaks a t tha t 
wavelength. This occurs first at optical and later at longer wavelengths. iGhisellini et al.l (|2007l ) 
ha ve found that the behavio r of 3 C 454.3 in 2005-2007 is consistent with the model suggested 
by iKatarzynski &: Ghisellini I (j2007l ). in which the dissipation site of an outburst depends on the 
bulk Lorentz factor and compactness of the perturbation propagating down the jet. Outbursts 
occurring closer to the black hole (BH) should have a more compact emitting region with a lower 
bulk Lorentz factor, F, and a stronger magnetic field, B. Greater compactness of the emission region 
intensifies the synchrotron flux as well as the high energy component produced via the synchrotron 
self-Compton (SSC) mechanism, while the external Compton (EC) high-energy component (inverse 
Compton radiation with seed photons from out side the jet) is suppress ed owing to a weaker Doppler 
factor, (5, resulting from the lower value of F. IGhisellini et al.l (120071 ) model the outburst in 2005 
as an event that occurred closer to the BH (F ~8, 6 ~13, B ~15 G, and size of the emitting 
region a ~5.5x 10' ^ pc) than the outburst in 2007 (F -16, 6 ^16, B ~9 G, and a ~8xl0"^ pc). 



Sikora et al.l (j2008l ) argue that the optical. X-ray, and millimeter light curves during the outburst 
in 2005 require a release of a significant fraction of the jet energy when the jet becomes transparent 
at millimeter wavelengths at the millimeter-wave "photosphere." These authors conclude that this 
phot osphere is located at ~10 pc from the BH, coinciding with the expected location of a torus of hot 
dust. ISikora et al.l (|2008l ) infer that the X-ray and 7-ray emission is most likely produced via the EC 
mechanism, with seed photons emitted by the hot dust scattered by relativistic electrons in a plasma 
with bulk Lorentz factor F ~20. Interpretations of 7-ray observations with AGILE" during autumn 
2007 and the densely sampled light curve provided by the Fermi Gamma-ray Space telescope 
starting in 2008 August, combined with simultaneous observations at longer wavelengths, involve 
higher-energy electrons that emit synchrotron radiation at near-IR and optical wa velengths as well 
as scatter external photon s from the broad line region to energies up to -^100 GeV (jVercellone et al 



2009 



Bonning et al.ll2009l ). 



Jorstad et al.l (|2005l ) monitored the quasar 3C 454.3 at 43 GHz with the Very Long Baseline 



Array (VLBA) bimonthly from March 1998 to April 2001, and determined parameters of the parsec- 
scale jet during this period based on the apparent speed of superluminal knots (also referred to as 
"components" of the jet) and time scale of their fiux variability. They found that the jet of 3C 454.3 
during this time span had physical parameters as follows: F=15.6ib2.2, Go=1.3°ibl.2°, (5=24.6ib4.5, 
and 0=O.8°ibO.2°, where Go is the viewing angle and 9 is the opening half-angle of the jet. The 



- 4 - 



Boston University group resumed VLB A monitoring of the quasar in June 2005, and has continued 
to monitor the source within a program of roughly monthly VLBA imaging of bright 7-ray blazars 
at 43 GHz. In this paper, we analyze disturbances seen in the quasar jet during the period of high 
optical activity in 2005-2008, and connect events in the jet with prominent variability at different 
wavebands. Observations cover the range of frequencies from 10^*^ GHz to 10^^ GHz, which provide 
significantly different angular resolutions at which an event is observed - 10-15 arcminutes at high 
energy frequencies, ~5 arcseconds at millimeter and sub-millimeter wavelengths, ~1 arcsecond in 
the optical and IR bands, and ~0.1 milliarcsecond (mas) with the VLBA at 43 GHz. This renders 
interpretation of multifrequency behavior challenging, and makes analysis of variability at different 
wavelengths along with VLBI images the main tool for understanding processes and mechanisms 
involved in the physics of blazars. 



2. Multifrequency Light Curves 

We have analyzed light curves of the quasar 3C 454.3 from 7-ray to radio wavelengths during 
the period 2004-2009. We use relative Juhan dates (RJD = JD - 2450000) to refer to the epochs 
of observation. 



2.1. Optical, Ultra- Violet, and near-Infrared Photometric data 



We use the optical light cu rve in R band collected by the WEBT collaboration ([Villata et al 



20061 . boOTlJRaiteri et ahlboOSal lbl) . which covers the period from 2003 December 28 to 2008 February 
8 (RJD: 3001-4505). To this we add the data from the 2 m robotic Liverpool telescope at the 
Observatorio del Roque de Los Muchachos (La Palma, Spain) and the 1.8 m Perkins telescope 
of Lowell Observatory (Flagstaff, AZ). We have extended the light curve up to 2009 January 
7 (RJD: 4839) using the data obtained at the Perkins telescope, 70-cm telescope of the Crimean 
Astrophysical Observatory (Nauchnij, Ukraine), 40 cm telescope of St. Petersburg State University 
(St. Petersburg, Russia), 2.2 m telescope of the Calar Alto Observatory (Spain), and 70 cm 
telescope of Abastumani Observatory (Republic of Georgia). At some telescopes, the photometric 
data have been carried out not only in R band but also in B,V, and / filters, although less frequently 
than in R band. We use these data for broad-band spectral analysis. The data from 2008 June 23 to 
December 10 (RJD: 4640-4810) are supplemented by measurements by the SMARTS consortium, 
posted at their websit^. 

We have performed the [/-band (3500 A) observations of the quasar with the Liverpool tele- 
scope from 2005 May 18 to Decen iber 10 (160 measur ements) and reduced the data in a manner 
used in the WEBT collaboration (iRaiteri et al.ll2008al ). We have obtained archived UVOT Swift 



""^ http://www.astro.yale.edu/smarts/glast/ 
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data from 2005 April 24 to 2009 January 27 in four filters, U (3501 A), UWl (2634 A), UM2 
(2231 A), and UW2 (2030 A) (86, 100, 89, and 95 measurements, respectively). The UVOT data 
were reduced using the HEASOFT 6.5 package and following the threads provided by the UVOT 
User's Guide and recommend ations contained in the release notes with parameters si milar to those 



used bylRaiteri et al.l (j2008al ). We have adopted galactic extinction values derived bv lRaiteri et al 
(|2008al ): 0.58 (?7), 0.73 {UWl), 1.07 (C/M2), and 1.02 mag {UW2). 



We use infrared (IR) JHK data collected at the 1.1 m telescope of the Main (Pulkovo) As- 
tronomical Observato ry of the Russian Academy of Sciences located at Campo Imperatore, Italy 
(jLarionov et al.ll2008l ). The values of Galactic extinction provided by the NASA Extragalactic 
Database are adopted to correct the photom etric estimates for Galactic absorption. We have 



applied the calibration of iMead et al.l (|l990l ) for all optical and IR measurements to transform 
magnitudes into flux densities. Figure [1] shows the optical light curve of the quasar in R band, the 
UV light curves in U and UWl bands, and the near-IR light curve in K band from January 2004 
to January 2009. 



2.2. Radio Light Curves 

The 230 GHz (1.3 mm) and 345 GHz (0.85 mm) light curves were obtained at the Submillimeter 
Array (SMA), Mauna Kea, Hawah from 2004 June 16 to 2009 January 31 (RJD: 3172-4862). 
3C 454.3 is a bright quasar included in an ongoing monitoring program at the SMA to determine 
the fluxes of compact extragalactic radio sources that can be used as calibrators at mm and sub-mm 
wavelengths. Data from this program are updated regularly a nd are availab l e at S MA websit^. 



Details of the observations and data reduction can be found in iGurwell et al.l (j2007l ) . 



The data at 230 GHz are supplemented by measurements carried out at the 30 m telescope 
of Instituto de Radio Astronomia Milimetrica (IRAM, Granada, Spain). The IRAM 30 m Tele- 
scope observed simultaneously at 86.24 GHz (3.5 mm) and 228.93 GHz (1.3 mm) by making use 
of the AlOO/BlOO and A230/B230 pairs of orthogonally linearly polarized heterodyne receivers, re- 
spectively. Every IRAM 30 m measurement was preceded by a cross-scan pointing of the telescope 
toward 3.5 mm and 1.3 mm calibration sources. Such measurements consisted of a series of wobbler- 
switching on-offs with total integration times of 4 min to 8 min, depending on the total flux density 
of the source and atmospheric conditions. Measurements of Mars and/or Uranus were obtained at 
least once per observing session in order to estimate and subtract residual instrumental polariza- 
tion, and to calibrate the absolute total flux density scale. Whenever these planets were not visible, 
the compact H II regions W3 OH, K3-50A, and NGC 7538, and/or the compact planetary nebula 
NGC 7027 were observed for calibration. The initial calibration of the amplitude was performed 



http://smal.sma.hawaii.edu/callist. html 
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through the telescope's onUne data processing procedures within the MIRA-GILDAS softwar^. 
The remaining data reduction involved the removal of outliers, an e levation dependent calibration, 
and the calibration of the absolute flux density scale as described in lAgudo et al.l (j2006l ) . The final 
flux density error of every measurement also included an additional 5% non-systematic uncertainty 
added in quadrature. Finally, the resulting data were averaged for those observing epochs on which 
more than one measurement was obtained. 

The 37 GHz (8 mm) observations were performed with the 13.7 m telescope at Metsahovi 
Radio Observatory of Helsinki University of Technology, Finland. The flux density calibration is 
based on observations of DR 21, with 3C 84 and 3C 274 used as secondary calibr ators. A detailed 
description of the data reduction and analysis is given in iTerasranta et al.l (|l998l ) . The radio light 
curves are plotted in Figure [H 



2.3. High-Energy Light Curves 



We obtained 63 measurements of 3C 454.3 with the Rossi X-ray Timing Explorer (RXTE) 
from 2005 May 11 to September 5 (RJD:3502-3619). Each RXTE visit lasted ~1 ksec. We used 
data taken with RXTE's Proportional Counter Array (PC A), which is collimated to have a FWHM 
1° field of view. However, there is a cat aclysmic variable star, IM Peg, located about 0.72° to the 
NW of 3C 454.3 (jPerrvman et al.l 119971 ). To completely eliminate contributions to the observed 
spectrum of 3C 454.3 from IM Peg, we chose an RXTE on-axis pointing po sition 0?52 to the S E 
of 3C 454.3; at this offset position, the PCA collimator efficiency is 45% (jJahoda et al.ll2006l ): 
we corrected for this when generating the PCA response matrices (see below). Reduction of the 
PCA data followed standard extraction and screening procedures, using HEASOFT version 6.0 
softwarfl. PCA STANDARD-2 data were collected from Proportional Counter Unit (PCU) 2 only. 
We use the "L7-240" background models, appropriate for faint sources. Response matrices were 
generated for each spectrum separately using the FTOOLS program PCARSP. We modeled the 
3-15 keV spectrum in XSPEC v. 11 with a single power -law and Galactic absorption corresponding 



to a hydrogen column density Ny{ = 7.2x10^^^ cm ^ ( Elvis et al.lll989l ). and calculated the light 



curve of X-ray flux density at 4 keV. 

For the period from 2005 May 10 to 2009 January 27 (RJD: 3501-4858) we used pubhcly 
available data from the Swift satellit^ that were carried out with the X-Ray Telescope at 0.3- 
10 keV (95 measurements) and processed through the XRT pipeline. We checked the data: the PC 
mode to WT mode transitions are well handeled, the rates are self-consistent, implying that there 
are no problems with the point-spread-function correction factors; the uncertainties are ^2%. 



^http://www.iram.fr/IRAMFR/GILDAS/doc/html/mira-html 
^http://heas arc. gsfc.nasa.gov/ftools/ 
■''http: / /www.swift. psu.edu/monitoring/ 



-7- 



We obtained Photon and Spacecraft data from the Fermi Science Support Center (FSSC) for 
observations with the Large Area Telescope (LAT) from 2008 August 5 to 2009 February 10 (RJD: 
4684-4873) within a 15° radius centered on the quasar. We then calculated the 7-ray light curve of 
3C 454.3 at 0.1-300 GeV with daily binning using the software and following the Analysis Threads 
provided by the FSSC. The procedure included selection of good data and time (programs: gtselect, 
gtmktime), construction of the exposure map for each day {gtltcuhe, gtexpmap), and modeling of 
the data by a single power law, optimized via a maximum-likelihood method (gtlike). We used 
the response function generated by the FSSC, and created a model file that consisted of a power- 
law model (prefactor and index) for 3C 454.3 and three other b right 7-ray sources in the field: 
PKS 2201-M71, CTA 102, and PKS 2325-^093 (|Abdo et al.ll2009al l. and diffuse emission models for 
Galactic and extragalactic background provided by the FSSC. We consider the quasar to be detected 
if the test-statistic, T5, calculated by the gtlike procedure exceeds 10, corresponding to at least 
a 3a detection level (jAbdo et al.ll2009al ). For such measurements we derived 7-ray flux estimates 
with uncertainties within 10-15%. We have s upplimented these 7-ray data by the measurements 
from RJD: 4649-4675 taken from Figure 1 of lAbdo et al.l (|2009bl ). which are not provided by the 
LAT Data Server. Figure [2] shows the high-energy light curves for three time intervals when the 
X-ray and 7-ray measurements were available along with the optical light curve. 



3. Polarization Observations 



We have obtained linear polarization measurements of 3C 454.3 at optical and mm wavelengths 
All o f these data were checked for consistency and corrected for statistical bias (jWardle &: Kronberg 
1973 ). 



3.1. Optical Polarization Data 



Our optical polarization monitoring in R band of 7-ray blazars, including 3C 454.3, began 
in 2005 May at the 70 cm telescope of the Crimean Astrophysical Observatory. The telescope is 
equipped with an ST-7 b ased photometer-pola rimeter. The details of the observation and data 
reduction can be found in lLarionov et al.l (|2008l ). 



Since 2005 September, we have obtained polarimetric observations at the Perkins telescope 
of Lowell Observatory (Flagstaff, AZ) with the PRISM camers[§ supplemented by a polarimeter 
with a rotating half-wave plate. These involve a series of 5-7 Stokes Q and U measurements for 
a given object. Each series consists of four measurements at instrumental position angles 0, 45, 
90, and 135° of the waveplate. Since the camera has a wide field of view (14' x 14'), we use field 
stars to perform both interstellar and instrumental polarization corrections. We use unpolarized 



^http: / /www.bu.edu / prism / 
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calibration stars from lSchmidt et al.l (|l992l ) to check the instrumental polarization, which is usually 
within 0.3%-0.5%, and polarized stars from the same paper to calibrate the polarization position 
angle. 

In 2006 November the 40 cm telescope of St. Petersburg State University joined the program. 
The telescope is equipped with a nearly identical photometer-polarimeter as that of the 70 cm 
telescope in Crimea, and the observations and data reduction are carried out in the same manner 
(jLarionov et al.ll2008l ). However, the polarization observations are performed without a filter, with 
central wavelength Agff ~ 670 nm. 

Startirig in 2007 July, the MAPCAT (Monitoring AGN with Polarimetry at the Calar Alto 
Telescopes ]I| program has made use of the Calar Alto Faint Object Spectrograph (CAFOS) in 
its imaging polarimetric mode on the 2.2 m Telescope at the Calar Alto Observatory (Almeria, 
Spain). Every i2-band polarization measurement consisted of four imaging exposures of ~50 s at 
0°, 22.5°, 45°, and 67.5° of t he A/2 plate in . CAFO S. Polarization measurements of at least two 
polarized standard stars from ISchmidt et al.l (| 19921 ) were performed to estimate the instrumental 
polarization, which is ~0.2%. The degree (p) and position angle (93) of l inear polarizati o n wer e 
obtained from their relation with the 8 flux measurements, as described in lZapatero et al.l (|2005l ). 



We have obtained polarization data at Steward Observatory with the 1.54 m Kuiper d uring 
a 2005 campaign from October 25 to November 4 (details given in iD'Arcangelo et al.ll2009l ). and 
within another, currently operating program at that telescope and the 2.3 m Bok Telescope to 
provide optical polarization measurements for bright 7-ray blazars from the Fermi LAT-monitored 
blazar lislo. The combined optical polarization data are presented in Figure [3j 



3.2. Millimeter- Wave Polarization Data 



The IRAM 30 m Telescope also performed linear and circular polarization at 86.24 GHz. The 
AIOO/BIOO pair of orthogonall y linearly polarized heterodyne receivers is connected to the XPOL 
polarimeter (iThum et al.ll2008l ). which provides 500 MHz bandwidth per receiver. We have carried 
out the initial calibration of the phase through the telescope's online data processing procedures 
within the MIRA- GIL DAS software. The subsequent polarimetric data reduction follows that of 
Thum et al.l (j2008l ) and lAgudo et al.l (j20ld ) . These procedures involve the removal of outliers and 
of systematic residual instrumental polarization still present in the data, calibration of the absolute 
flux density scale, estimation of remaining non-systematic or statistical uncertainties still affecting 
the data, and quadratic addition of the latter to the final uncertainties of the measurements. Such 
final errors are ^ 0.5% in Stokes Q and U parameters, and ~ 0.3% in Stokes V parameter. The 
data are plotted in Figure [3l 



http://www.iaa.es/~iagudo/research/MAPCAT/MAPCAT.html 
*http://james. as. arizona.edu/~psmith/Fermi 
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4. VLBA Observations 

We observed 3C 454.3 in the course of a program of monthly monitoring of bright 7-ray blazars 
with the VLBA at 43 GHz (7 mmj^. During the period from June 2005 to May 2009, we obtained 
44 total and polarized intensity images of the quasar at a resolution of ~0.3x0.1 milliarcseconds 
(mas). Figure [4] presents a sequence of total intensity images, with several images per year convolved 
with the same beam and with co ntours based o n the same global map peak. We performed the 
data reduction in the manner of 



Jorstad et al. 



(200 



sl) using the Astronomical Image Processing 



System (AIPS) and Difmap (IShepherd I Il997l ). In addition to the standard resolution beam of 
uniform weighting of the uf-data, we have convolved the images with a beam of FWHM diameter 
of 0.1 mas comparable to the resolution of the longest baselines of the VLBA. Figure [5] shows the 
morphology of the parsec-scale jet of 3C 454.3, convolved with a standard beam along with an image 
with the smaller beam that features knots Kl, K2, K3, and C found in the jet (see §5). Although 
the synthesized beam of the VLBA observations along the jet is ~0.1 mas, we have convolved the 
core region of the jet with a circular beam with a FWHM diameter of 0.05 mas for epochs when all 
antennas were in operation and the weather was favorable. Since this is only 30% smaller than the 
FWHM angular resolution of the longest baseline, 0.07 mas, the array partially resolves structures 
on 0.05 mas scales. At the majority of suitable epochs, the resulting super-resolved map of the 
core region consists of three components, as shown in Figure [H Such a structure is especially 
apparent when the core is not disturbed by passage of a moving knot [bottom panel). To support 
the conclusion that this fine structure in the core is real, we have constructed similar super-resolved 
images of the quasar CTA102, obtained with the VLBA along with 3C 454.3. Since CTA102 is 
close to 3C 454.3 in the sky, the observations have similar uf-coverage. The resulting images of 
CTA102 contain only single-component cores. This suggests that the fine structure of the core of 
3C 454.3 is not an artifact of the smaller beam. 

Since VLBA data yield only the relative orientation of the electric vector position angle 
(EVPA) across a source, special attention was paid to the calibration of the absolute EVPA values. 
We employed data from the NRAO VLA/VLBA Polarization Calibration Pag which provides 
EVPA measurements with the VLA at 43 GHz, integrated over the entire source for several ob- 
jects from our sample observed with the VLBA (OJ287, 1156+295, 3C 273, 3C 279, BL Lac, 
and 3C 454.3). For other epochs witho u t con temporaneous VLA observations, we used the D- 



terms method discribed bv iGomez et al.l (|2002l ). The accuracy of the absolute EVPA calibration 
is within 4-8°. iJorstad et al.l (|2007l ) provide an estimate of the Faraday rotation measure [RM) of 
(-6.8 ± 3.7) X 10^ rad m ^ in the mm-VLBI core of 3C 454.3, corresponding to a correction for 
the EVPA at 43 GHz of -19° ± 10°. Although the estimate of RM has a large uncertainty and 
is not simultaneous with the VLBA observations reported in this paper, we find that applying the 
RM correction at 43 GHz improves the similarity between the behavior of EVPAs at 86 GHz and 



http://www.bu.edu/blazars/VLBAproject.html 
^''http: / /www. vla.nrao.edu / astro / calib/polar / 
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43 GHz. Figure [3] shows the results of our measurements of degree (^43) and position angle (9543) 
of polarization in the core region of the quasar. 



4.1. Modeling of Images 

We have modeled our calibrated uu-data with point source brightness distributions using the 
Difmap task MODELFIT. The number of point-like components required to fit the data was de- 
termined by the best agreement between the model and data according to values, with reduced 
ranging from 1.0 to 5. Since we have roughly monthly observations with the VLBA, the model 
yielding the best agreement at a given epoch was used as an input model for the following epoch 
in order to maximize consistency across epochs given the complexity of the brightness distribu- 
tion. We then performed 100 iterations to adjust the input model to the wf-data of the second 
epoch, editing the model if deleting or adding components was required to represent the image. 
We assumed that the core is a stationary component at relative right ascension and declination 
(0,0), and then determined for each component the total and polarized flux intensity, S43 and S"^^, 
respectively, distance from the core, R, position angle relative to the core, 0, degree of polarization, 
P43, and position an gle of polarizat i on, (^4 3. We have identified components across the epochs in the 



manner described in lJorstad et al.l (|2005l ) and used their method, as well as the same cosmological 
parameters {Q^n = 0.3, Q\ = 0.7, and Ho=70 km s~^ Mpc~^), to calculate the proper motions 
and apparent speeds of moving features. Figure [1] shows the light curve of the VLBI core, while 
Figure [7] presents the positions of all components (except the core) brighter than 50 mJy found 
within 1 mas of the core during the period 2005 June to 2009 May (RJD: 3547-4982). 

Our analysis of the VLBA images concentrates on the innermost region of the jet. Figure [8] 
shows the projected jet direction, Gjet, determined at each epoch according to the position angle 
of the brightest component within 0.1-0.3 mas of the core. The average jet direction is {Qjet)=- 
95°ib8°. Figure [8] shows that there are significant deviations of 0jet from the average, which might 
reflect intrinsic changes of the jet direction or shifts in the brightness distribution across the width 
of the jet. 



Disturbances in the Inner Jet 



Analysis of the images reveals three features moving with respect to the core for which the 
apparent speed, /3app, and time of ejection (coincidence with the core). To, can be estimated. We 
have designated the features as components K1,K2, and K3 (Fig. [7]). The majority of images 
show a brightness enhan cement at 0.6-0.7 ma s from the core, which is most likely associated with 
component C reported by lJorstad et al.l (j2005l ) as a stationary knot at position R ~ 0.63 mas. Knots 
Kl and K2 move ballistically only within 0.2 mas of the core. Beyond 0.2 mas, Kl accelerates 
and fades, while K2 decelerates and brightens. Figure [9] shows light curves and trajectories of 
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components. Although the separation of K3 from the core is at the Umit of our resolution, modeling 
requires the presence of two components in the core region at all epochs starting from 2008 June, 
with a persistent increase in separation between the components at later epochs. Figures \W\ \TT\ 
and [12] show sequences of total and polarized intensity images displaying the jet evolution that we 
attribute to the emergence of components Kl, K2, and K3, respectively. Knots Kl and K2 are 
distinct on the high-resolution total intensity maps and have better alignment of the position angle 
of the polarization with the jet direction than does the polarization of the core. The brightest knot, 
K3, fades dramatically while it is still very close to the core. 

Since we are interested mainly in the time of ejection of the components and their speed in the 
vicinity of the core, we have used only those epochs when the components moved ballistically (the 
first 10 and 8 epochs for Kl and K2, respectively, and all 11 epochs for K3) for deriving kinematic 
parameters. The p arameters and their uncertainties were obtained in the manner described in 



Jorstad et al.l (120051 ') 



We modeled components by point sources because such an approach gives better consistency 
between epochs than, for example, modeling by components with circular Gaussian brightness 
distributions. This suggests that the knots are very compact, with angular size smaller than half 
the beam axis along the jet, ^0.07 mas. We have adopted a FWHM size of 0.05 mas for all three 
components to estimate Doppler and Lorentz factors and viewing angle indep endently by using 



the method developed for high-frequency VLBI monitoring (I Jorstad et al.ll2005l ). Table [T] lists for 
each knot the time of ejection, To, the highest flux density measured, S^^, the epoch of 5™^^, 
proper motion, /i, apparent speed, /3appi Lorentz, F, and Doppler, S, factors, and viewing angle, Go- 
There is a significant difference in the viewing angles of components K2 and K3, and the Lorentz 
and Doppler factors of K3 are twice those Kl and K2, as well as twice th ose derived from th e 



parameters of several components detected in 3C 454.3 between 1998-2001 (I Jorstad et al.ll2005l ). 



Figure [8] shows that the most dramatic change in the projected jet direction o ccurred after e j ection 



of component K3. If we assume that the average jet parameters derived by iJorstad et al.l (|2005l ) 
give the most probable values of viewing and opening angles (Go=1.3° and 6=0.8°, respectively) 
then component K3 with 0o=O.2° should occupy the side of the jet closest to the line of sight, while 
component K2 with 0o=2.5° should move along the far side of the jet. Figure [8] supports such an 
interpretation, since, after the ejection of component K2, the projected jet direction swings to the 
north, while after the ejection of K3, Qjet turns to the south. Component Kl appears to follow a 
path closer to the jet axis. In addition. Figure [8] shows a sharp swing of Qjet near RJD~3750-3800, 
possibly indicating the emergence of a new knot from the core that we cannot distinguish from Kl, 
which is bright and close to the core at this time (see Fig. [7]). 
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6. Correlation Analysis between High Energy Variations and Light Curves at 

Optical and mm Wavelengths 

The X-ray observations were carried out less uniformly than observations at other wavelengths. 
There are three intervals when X-ray observations were obtained: I - RJD: 3502-3619, II - RJD: 
4292-4450, and III - RJD: 4613-4858 (see Fig. [2]). We have performed correlation analysis be- 
tween the X-ray, optical, and 230 GHz light curves separately for these intervals because (i) there 
are significant gaps in X-ray coverage between the intervals, (ii) X-ray measurements at different 
intervals are obtained with different instrume nts, and (iii) correlatio n between X-ray and longer 
wavelength light curves can change with time (IChatterjee et al.l 120081 ) . For interval III, we include 
the Fermi 7-ray light curve in the correlation analysis. Al l light curves are group ed into 1-day bins. 
We employ the code developed for correlati on analysis (IChatteriee et al.ll2008l ). which calculates 
the discrete cross-correlation function, CCF ([Edelson &: Krolik Ill988j). finds the peak of th e CCF, 
fmax, and determines the position of the centroid near the peak (| White &: Peterson 1 1 19941 ). This 
yields an estimate of the significance of the peak as well as the delay, r, between variations at 
different wavelengths if the peak is significant at a confidence level 0.05. 

Interval I has the best X-ray data sampling, provided by RXTE monitoring. We have correlated 
this X-ray light curve with the optical and 1 mm light curves over exactly the same period, while 
performing the correlation between the R-band and 230 GHz light curves over a longer period, 
RJD: 3499-3750, which has similar sampling at the two wavelengths. Figure [13] {left panel) gives 
the results of the CCF analysis, revealing very strong correlation between variations at X-ray 
and longer wavelengths. During interval II, there are only 21 measurements obtained with Swift, 
insufficient for correlation analysis. We compute the CCF only between the optical and 230 GHz 
light curves for the period RJD: 4200-4505. Figure [13] (right panel) shows a moderate correlation 
between variations at the two wavelengths. During interval III, the X-ray (Swift), optical, and 
230 GHz light curves have good coverage, with 72, 137, and 115 points, respectively; daily 7-ray 
fiuxes became available starting at RJD: 4684. We hav e added 7-ray measurements from RJD: 
4649-4675 presented in Figure 1 of ]Abdo et al.l ()2009bl ) to the 7-ray light curve for correlation 
analysis between 7-ray and optical variations. In the case of the X/7-ray correlation, the 7-ray 
data during RJD: 4649-4675 are not used due to the absence of X-ray measurements during this 
period (see Fig. [2]). Figure [T^ {left panel) presents the cross-correlation functions between X-ray, 
optical, and 230 GHz light curves, while the right panel shows the CCF between the 7-rays and 
lower frequencies. Table [2] summarizes the results of the correlation analysis for all intervals and 
wavelengths. 

Table [2] shows a strong correlation between optical and 230 GHz variations with a delay of 
mm-wave emission with respect to the optical variations. Time lags between the optical and 1 mm 
light curves decrease significantly, from ~50 days to ~10 days, from interval I to III, respectively. 
Raiteri et al.] ()2008a] ) have reported a change in delay between the optical and 230 GHz light curves 
of the quasar during 2007-08 relative to 2005-07, consistent with the difference in ropt/230 that we 
find between intervals I and II. During interval III, the delay became even shorter. Analysis of 
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the X-ray/optical correlation suggests that variations at the two wavelengths are simultaneous 
within 1 day. However, a change from positive (optical variations lead) to negative (X-rays lead) 
delay might have taken place from interval I to interval III, although we judge this to be of low 
significance given the uncertainties of the delays. Nevertheless, such a possibility is supported by 
the X-ray/230 GHz delays as well: during interval I, the X-ray/230 GHz delay is shorter than the 
delay between the optical and 230 GHz light curves; the opposite situation (|Topt/23ol > kx/23ol) is 
observed for interval HI, but again the differences are within the uncertainties of r. 

Interval HI is especially interesting, since it contains the Fermi 7-ray light curve. The 7-ray 
variations show the best correlation with the optical ligh t curve, with no dela y exceeding 1 day 



(Fig. [HI right panel), in agreement with the findings of iBonning et al.l (|2009l ). However, these 
authors reported the absence of a correlation between 7- and X-rays, while our analysis indicates 
the presence of such a correlation, although the peak of the CCF is less prominent than that between 
the 7-ray and optical light curves. From this analysis, we conclude that the 7/X-ray variations 
are simultaneous within 2 days. Our X- and 7-ray light curves include additional data after RJD: 
4750, which probably explains the discrepancy between our correlation result and that obtained by 



Bonning et all ((20091). 



We have performed correlation analysis between the 230 GHz and 37 GHz light curves using all 
data from RJD: 3000 to RJD: 4850 (see Fig. [T]). The CCF between the two wavelengths presented 
in Figure [15] shows a global maximum with a delay r23o/37=— 215ib30 days and two local maxima 
at r23o/37=— 73±30 days and T23o/37=0±10 days. This is connected with a change of the delay 
between the variations at the two wavelengths with time. Indeed, the CCF for the period from 
RJD: 4000 to RJD: 4850, which excludes the largest mm-wave outburst, gives a well defined peak 
at T23o/37=0=b7 days (Fig. [T5l dotted line). The change of the delay between the 230 GHz and 
37 GHz light curves from ~200 days to days reflects the transition of mm-wave emission of the 
source from optically thick to optically thin (see §8), which also leads to a decrease of the delay 
between optical and mm-wave variations. 



7. Comparison of Timing of Flux Variability and Ejection of Superluminal 

Components 

We decompose the light curves at optical and r adio wavelengths into individual flares with expo- 
nenti al rise and decay by applying the method of Valtaoja et al. (j 19991 ) (see also Chatter jee et al, 



20081 ). The i?-band light curve is smoothed with 1-day binning, while the mm-wave curves are 
smoothed with 5-day binning. The method gives the epoch of peak flux, T^, output power, Ei, 
(area under the curve fitting a flare), and width, Wi, (average of the rise and decay times). The 
decomposition yields 10 outbursts in the i?-band light curve and five outbursts in the mm-wave 
curves, each with >10 days (Table [3|). We consider these to be the most reliable flares, based on 
the sampling of the data. We have cross-identified outbursts at different wavebands by assuming 
that events associated with each other should have a difference in times of the peaks correspond- 
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ing to the delay found from the correlation analysis of light curves at the two wavelenths, i.e. 
Tj^^ — Ti,2 ^ Tui,u2 (see §6). We also analyzed relative values of the parameters Ei, and Wu to find 
additional support for the cross-identification. This implies, for example, comparison of Eopt/E^^^ 
with £'23o/-£'230^ ) where E^^ and -E'230^ maximum output powers derived among optical 

and mm- wave outbursts, respectively. The results of cross-identification are shown in Table [3l We 
have associated the optical outbursts peaked at RJD: 3501, RJD: 4301, and RJD: 4437 with the 
230 GHz outbursts with maxima at RJD: 3534, RJD: 4330, and RJD: 4454, respectively, because 
differencies in the peaks (—33, —29, and —17 days) agree very well with the delays found between 
the two wavelengths (Table [2]) . The rest of the optical flares do not have counterparts at mm- 
wavelengths. All 230 GHz outbursts have corresponding 37 GHz outbursts, based on the results of 
the correlation analysis. 

We have also related the mm-wave outbursts to the ejection of superluminal knots, based on 
the assumption that an outburst and an ejection are connected if the latter occurs within the 
width of the outburst. Table S] gives the delays between the time of the ejection and the peak of 
the corresponding outburst at different wavelengths, ATi, = Tu — Tq. It shows that each ejection 
occurred within 20—30 days of the peak of a 230 GHz outburst. We denote the 230 GHz outbursts 
as Mki, Mk2, and Mks in accordance with the designation of related superluminal component 
(Fig. [1]). We then associate a newly appearing superluminal knot with an optical outburst if the 
latter is related to a mm-wave outburst. As a result, the three major optical outbursts with the 
highest power are associated with newly emerging superluminal components. We denote these 
optical outbursts as Oki, Ok2, and Ok3 in a similar manner as the mm-wave outbursts connected 
with the ejections (see Fig. [T] and Table [3]). 

Analysis of Tables [3] and S] reveals that (i) all ejections trigger an increase in optical and 
mm-wave activity; (ii) the connection between optical outbursts and ejections of superluminal 
components is not one to one, since either one ejection can trigger multiple optical oubursts or 
there are optical outbursts that are not connected with a new knot; an d (iii) all ejection s occu r 



during the rising branch of a 37 GHz outburst, as previously reported by ISavolainen et al.l (120021 ). 
In the case of the Kl event there is a pattern that the strongest optical outburst is followed by the 
230 GHz outburst and by the ejection of a superluminal knot from the 43 GHz core. The delays 
between the peaks of the outbursts and time of ejection of the knot can be interpreted in terms 
of spatial separation of sites of the outbursts with respect to the core. In the case of events K2 
and K3, such an interpretation is possible as well, taking into consideration uncertainties in AT^. 
However, another interpretation, which better corresponds to the derived delays and reflects the 
differences in properties of Ok2 and Oks with respect to Oki, suggests that the sites of both the 
optical and 230 GHz outbursts are located in the mm-wave core and that the outbursts occurred 
during the ejection of a superluminal knot. A delay between the peaks of associated optical and 
230 GHz outbursts can be caused by the difference in time when the maximum number of relativistic 
electrons radiating at optical and mm- waves was achieved: the time of the maximum number of 
relativistic electrons accelerated to optical-emitting energies coincides with the ejection time, while 
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acceleration of electrons to mm-wave-emitting energies continues for 20-30 days longer. 



The extensive optical data collected by the WEBT collaboration (jRaiteri et al.ll2008al ) contain 
a number of epochs with several measurements during a single date. These authors have noted 
strong intraday variability in 3C 454.3 in the optical R band on a time scale of tvar ~l-5 hr in 
December 2007. We have calculated the means and standard deviations for R magnitudes for 
each Julian date with 7 or more measurements. Figure [16] presents the standard deviations, 
vs. epoch, along with the average uncertainty for individual measurements on a given date, ctR 
. There is possible intraday variability during the ejection of component Kl, and pronounced 
intraday variability during the emergence of knot K3. According to our estimations, component 
K3 has the highest Doppler factor that we observed in the jet of 3C 454.3 from 2005 to 2009. 
This supports the idea that the knot is responsible for the short time scale of the variability 
seen at optical wavelengths. The angular size of the region that can prodice such a variability, 
Ovar ^ cAtvar(l + z)/{5Dp^) f» 10~^ arcscconds, where Da is the angular distance to 3C 454.3. 
The size of that we adopt from th e VLBA imaging is ~10^ larger than a^r- H owever, analysis 



of polarization properties in quasars (jJorstad et al. 



2007 



D'Arcangelo et al.l 120071 ) shows that the 



polarized optical emission, although partially co-spatial with the mm-wave core, occupies a volume 
that is thousands of times less than does the polarized emission at 43 GHz. Therefore, it is possible 
that the site of an optical flare lies within a superluminal knot passing through the core, and is 
confined to such a small region in the knot that it can vary on extremely short timescales. 



8. Spectral Behavior 

We use 340 GHz, 230 GHz, 86 GHz, and 37 GHz data to calculate mm-wave spectral indices, 
ttmm {Sy oc u^'^). The spectral index for a given Julian date was calculated if there were mea- 
surements at a minimum of three wavelengths within 2 days of the date. When there were several 
measurements over a short time period, the observations nearest to each other were used. Fig- 
ure [T7] reveals that the mm-wave spectral index exhibited significant variability. The uncertainty 
in amm reflects the scatter about the best-fit power law. During a quiescent state amm is typical 
of a compact, flat-spectrum radio source, amm =0.18±0.04, where Omm denotes a quiescent state 
chosen during the period RJD: 3900-4200 (Fig. [1]) and denoted by QS. The spectrum became 
strongly inverted during outburst Oki and slightly inverted during outburst Ok2- The spectral 
index reaches a local minimum at epochs of superluminal ejections. After an ejection, Omm changes 
sharply, rising to the value of Omm, and then continuing to increase up to a value amm = 0.5-0.6, 
signifying that the emission is close to being optically thin. The epochs of the two global maxima 
of 

Q^mm coincide with the peaks of the sharp mm-wave outbursts at RJD: 3781 (-A'^l) and RJD: 
4705 (M2), which are identified neither with optical outbursts nor with ejections of superluminal 
components (Table O Fig[T]). The peaks of Ml and M2 do, however, coincide with the epochs of 
maximum 43 GHz flux of knots K\ and i^3, respectively (Table [1]). 
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Raiteri et all (|20oi boOSal lbll plot the spectral energy distributions (SEDs) of the quasar at 



different brightness levels. These authors interpret the behavior of the optical- UV part of the spec- 
trum as a result of superposition of beamed jet emission, blended line emission from the broad 
line region (little blue bui np), and emi ssion from the accretion disk (big blue bump), which agrees 



with previous findings by ISmith et al. I (11988). The compound SED explains the flattening of the 



spectrum with a decrease in brightness, since at a low emission state the contribution of t he blue 



bumps becomes more pronounced. We use the method suggested by iHagen-Thorn I (|l997l ) to de- 
termine the spectral characteristics of a component of the source responsible for flux variability 
observed on timescales of days and weeks. Such a component most likely originates via the syn- 
chrotron mechanism and has a power-law spectral energy distribution from IR to UV wavelengths. 
The method allows one to determine spectral indices of the synchrotron emission separately from 
the thermal emission, which varies on much longer timescales. The detailed description of the 
m ethod and its applic a tion t o multifrequency optical and near-IR data for the event Oki are given 



m 



Hagen-Thorn et al.l (|2009l ). We have performed a similar analysis for other events. Table [5] gives 
the near-IR and optical spectral indices of synchrotron components during outbursts Oki, Ok2 + 
Ok^-i and M2, as well as the quiescent period QS. Outbursts Ok2 and Ok^ are combined owing 
to the small number of multifrequency optical and IR observations obtained separately for each 
outburst. 

We have derived spectral indices of synchrotron components in the UV-region, applying the 
same method to multifrequency Swift UVOT measurements. We have constructed relations between 
fluxes in {/, UM2, and UW2 bands relative to fluxes in UWl band, with the largest number of 
measurements, during two periods RJD: 4292-4450 {Ok2 + Okz) and RJD: 4613-4858 (M2), when 
multifrequency UV-data were available (see Fig. [T]). Figure [TSl shows the flux-fiux dependences, 
which can be fit by straight lines, Si = A\ + i?i5uwi- The slopes of the regressions, are used to 
construct the relative SED of the synchrotron component responsible for the variability. Figure [19] 
shows the relative SEDs of synchrotron components at the UV bands along with the those obtained 
for IR and optical parts of the electromagnetic spectrum (Table [5]). The relative UV SEDs are very 
similar during both events. However, they cannot be represented by a single power-law owing to 
significant flattening of the spectrum at wavelengths below 2630 A. We have determined the spectral 
indices for the soft UV (2634-3500 A), a^^^, and for the hard UV (2030-2634 A), a^^^, emission 
separately (Table [5]) . We stress that this hardening of the UV spectrum cannot be connected to the 
contribution of the blue bumps to the SED owing the short timescale of variability of the emission. 

Figure [19] and Table [5] show that the optical spectra of the variable components are similar at 
the flaring st ates and flatter than a^p^^ . This is a characteristic of optical variability often observed in 



blazars (e.g.. lLarionov et al.ll2008l ): variable sources are bluer at brighter flux levels. In the quiescent 
state a^j^ is flatter than a^p^. by ~0.5, as expected for a synchrotron source with relativistic electron 
energy distribution N{E) = NoE~^'^ when optically emitting electrons suffer strong radiative losses. 
The signiflcantly steeper with respect to a^J^ during event Oki is unusual. This effect can 
be caused by undersampling the IR-data during the brightest stage of optical emission (see Fig [T]) . 
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However, if it is real, a hardening of the synchrotr on spectrum can occur at electron energie s 
where Klein-Nishina (KN) effects become important (jPermer &: Atovan II2002I : iModerski al.ll2005l ). 
if inverse Compton losses dominate over synchrotron losses, i.e., Uphot ^ -^^/(Svr), where Uphot 
is energy density of the photon field and B is strength of the magnetic field. Such a condition 
might occur during the dramatic outburst Oki that yields diminishing inverse Compton losses at 
frequencies u > 2.6 x 10^^6B/{1 + z){v'' /lO^^)^'^ Hz, where z/' ~ Tuseed and Vseed is the frequency of 
the external photon field. If we assume that 6 and T correspond to knot Kl and I'sced ~ 5 x 10^^ Hz, 
a hardening of the synchrotron spectrum at optical wavelengths, v > A x 10^^ Hz, requires the 
magnetic field B ~4.5 mG, while for B ~1 G the external photon field should peak in U-band, 

I/seed ~ 7.5 X 1014 Hz. 

The soft UV spectrum is stee p, with g^Jy s a2.3, an even steeper synchrotron spectral index 
at 2500 A, ~ 2.7, was derived by ISmith et al.l (|l988l ) from multifrequency polarization data of 
3C 454.3. The hardening of the synchrotron UV spectrum beyond 2500 A that we observe during 
events Ok2+Ok3 and M2 can be caused by inverse Co mpton losse s in the Klein-Nishina regime, 
as described above. According to the SED presented by lAbdo et al.l ()2009bl ). the inverse Compton 



peak exceeds the synchrotron peak during M2 event by at least a factor of ^^5, w hich implies 



that the conditions for inverse Compton losses should be achieved (jModerski al.l 120051 ). For events 
Ok2+O k?, and M2 the .jet p arameters are uncertain. We adopt the average jet parameters r=16 and 
5 = 25 (jjorstad et al.ll2005l ) that result in a hardening of the synchrotron spectrum at > 2 x 10^^ 
if z^seed ~ 3 X 10^^ Hz and B ?»50 mG or \i B G and i/seed ~ 2 x 10^"^ Hz. 

X-ray spectral indices at 2.4-10 keV were calculated in the process of the RXTE data reduction. 
We use the X-ray spectral index from the epoch closest to the time of the peak of Oki to characterize 
the outburst. For outbursts 0x2, 0x3, and M2 we have derived at 2.4-10 keV from Swift data 
obtained at epochs closest to the dates of the peaks of the outbursts using the same spectra l mode l 
as for the RXTE data,. The X-ray spectral index for the quiescent state is from lRaiteri et al.l (|2007l ). 
derived from XMM-Newton measurements on 2006 December 18/19. Table [5] summarizes the X-ray 
spectral indices and shows that during outbursts Ok2 and Oks Ox is similar to that in the quiescent 
state, while during the most powerful optical outburst, Oki, Ox became steeper than a^. 

We have determined 7-ray spectral indices, ao.i-300Get;i using a single power-law model. Fig- 
ure [20] shows variations of the 7-ray spectral index during RJD: 4684-4873. The average value 
(ao.i-300Gev)=l-46±0.29 corresponds to the average o f the optica, l and ne ar-IR spectral indices 
of the synchrotron component during outburst M2. lAbdo et al.l (j2009bl ) have found that the 
7-ray spectrum of 3C 454.3 during M2 event is described by a broken power-law model with 



a. 



low_ 



-1.27±0.12 and aij'^ =2.50±0.35 and a break at 2.4±0.6 GeV. This implies a close correspon- 
dence between al^™ and and between ctij'^^ and the spectral index of soft synchrotron UV 



emission, ctyy^ (see Table [5]). 
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9. Polarization Behavior 

Figure [3] shows the parameters of the polarization from the whole source at both optical 
wavelengths and 86 GHz, and at 43 GHz in the VLBI core region. The data reveal a very 
wide range of variability of the degree of optical polarization, from ~0% to 30%. The range and 



timescale of the variability decrease with wavelength, in agreement with the findings of lJorstad et al 



(|2007l ). The average values of the degree of polarization are (popt)=7.4ib5.3%, (p86)=2.3±1.4%, 
and (p43)=1.6=b0.8%. At all three wavelengths the position angle of polarization covers the entire 
range from —180° to 0°, although (popt varies more rapidly than that of the mm emission. We 
have collected 253 simultaneous measurements of the optical flux and degree of polarization, which 
produce a statistically significant correlation between these parameters (linear coefficient of correla- 
tion, />=0.32). The connection between the fiux and degree of polarization is weaker in the 43 GHz 
core (p=0.21) and absent at 86 GHz (p=0.05). 

We have constructed distributions of alignment of EVPAs at optical and 86 GHz wavelengths 
from the whole source, and at 43 GHz from the VLBI core region with respect to the jet direction. 
The distributions are presented in Figure [2T] for the all measurements at each frequency, as well as 
for cases when the degree of polarization is higher than the average at a given frequency. The optical 
and 43 GHz core EVPAs extend across all possible directions with respect to the jet axis, although 
(/7opt tends to be perpendicular to the jet, while ip^^ aligns more often with the jet direction. These 
tendencies are more prominent when the degree of polarization is high. This might be partially 
caused by the larger uncertainties in EVPA when p is low. The EVPA at 86 GHz maintains good 
alignment with the jet direction independent of degree of polarization. However, this evidence of 
better alignment may be biased by the poorer time coverage at 86 GHz since, in general, Figure [3] 
shows good agreement between the EVPA measurements at 86 and 43 GHz. 

Table [6] presents the optical polarization position angle and EVPA in the 43 GHz core for es- 
sentially simultaneous observations (within a week, less than the time scale of significant variability 
at 43 GHz). The optical and 86 GHz polarization position angles measured within 2 days of each 
other are listed in Table [71 There are 19 simultaneous pairs of optical and 43 GHz measurements, 
with 8 pairs (42%) showing alignment between (/9opt and (^943 within the Icr uncertainty of the mea- 
surements. The optical and 86 GHz EVPAs agree in only 2 out of 12 pairs (17%). Alignment 
between (/?opt and ip^s is observed independently of direction of polarization with respect to the jet 
axis, while two cases of agreement between ipopt and pse occur when the EVPAs are close to the 
jet direction. If optical and 43 GHz EVPAs are independent, Bayes theorem yields a probability 
of <3xl0^^ that the number of observed alignments between the optical and 43 GHz polarization 
position angles occur by chance. This implies that the polarized emission at different wavelengths 
is govern ed by common properties or processes ra. t her than by chance, in a greement with previous 



findings (jGabuzda et al.1 12OO6I : Ijorstad et al.l l2007l : iD'Arcaneelo et al.l 120071 ') . 
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9.1. Rotation of Position Angle of Optical Polarization 



Analysis of the optical polarization behavior during optical outbursts shows that each ejection 
coincides with a rotation of the optical polarization position angle (Fig. [22]) . The rotations have 
similar rates for events Kl and K2, v^ = 8.7±1.1 and v^— 7.7^:0.8 degrees per day, while the 



rotation during the ejection of K3 is slower, z;^t=5.2ib0.7 degrees per day. We interpret this as 
evidence that the rotation of ipopt originates in the acceleration and collimation z one (ACZ). We 



can rule ou t a stochastic rotation of the position angle due to plasma turbulence (jMarscher et al 
20081 . 12OIOI ) , given that all three events have the same sense of rotation and coincide with optical 



flares. In the ACZ, Poynting flux is converted to bulk kinetic energy through the Lorentz force, 
which r otates the MHD outflows a nd forms helical trajectories of disturban ces propagating down 
the jet (jviahakis fc Konigl I I2OO4I : Iviahakis I bood : Komissarov et al.1 boOTI ). The Lorentz factor 
should increase linearly with cross-se ctional radius r , while the speed of rotation should decrease as 
to conserve angular momentum ( Vlahakis II2006I ). If the polarized optical emission is produced 
in the ACZ by an emission feature that follows streamlines that execute a spiral path about the 
axis, then the speed of rotation of different disturbances should be inversely proportional to the 
final bulk Lorentz factor. The jet flow must travel farther from the BH to attain a higher Lorentz 
factor, hence the size scales of the site of EVPA rotation should be larger and the position of the 
core should shift downstream relative to states with a lower terminal Lorentz factor. In support 
of these predictions, the lowest rate of rotation is observed during event K3, which corresponds to 
the highest value of F, while v^l ~ v^, in agreement with Fki ~ Fk2 (see Tabled]). Furthermore, 
there is a discontinuity in the motion of knot K2 near the time of ejection of K3 (see Fig. [7]), which 
could be explained if the 43 GHz core shifted downstream at this time. At the end of the optical 
EVPA rotation in events K2 and K3, the EVPA became chaotic, with average value perpendicular 
to the jet (see Fig. [22]) . This is consistent with an optical outburst arising from an interaction 
between a disturbance propa gating through a turbul ent section in the jet and a standing shock 
system in the mm- wave core (jPArcangelo et al.ll2007l ). 



Rotation of the optical EVPA associ ated with a disturban ce propagating down the jet has 



been observed p reviously in BL Lacertae (iMarscher et al 



PKS 1510-089 (Larionov et al. 



2008 



Marscher et al 



20081 ) and in the quasars 3C 279 and 



2OIOI ). There are several other properties of 



such rotations that are similar in BL Lac, 3C 279, PKS 1510—089, and 3C 454.3: (i) at the beginning 
of the rotation, the degree of polarization is low (~l-2%); (ii) maximum fractional polarization is 
achieved when the E-vector aligns with the jet direction (for event Kl, Pinax=22%, V'max=— 77°; 
for K2, Pjnax=26%, V'max=— 87°; event K3 is undersampled) ; and (iii) at the end of the rotation, 
the degree of polarization decreases. We infer that a rotation of the optical polarization position 
angle during a prominent optical outburst is a common feature of blazars. In addition, in 3C 279 
both the optical and 43 GHz core EVPAs rotate as a superluminal knot emerges from the core 
(jLarionov et al.ll2008l ). A few measurements of the EVPA in the core of 3C 454.3, simultaneous 
with optical polarization observations, agree with the optical EVPAs (Fig. [22]) . However, the 
agreement occurs when ipopt and 9343 are closely aligned with the jet direction (Table [6]), so that it 
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is possible that the core EVPA remains near that preferred angle instead of rotating with (popt- 



Despite agreement between ipopt and (^43 for a number of essentially simultaneous observations 
(see Table [6]), the majority of the high optical polarization measurements occur when the optical 
EVPA is perpendicular to the jet, while at 86 GHz and in the 43 GHz core, the majority of the E- 
vectors at high degrees of polarization lie along the jet direction. There are multiple possible reasons 
for this discrepancy: (i) high mm-wave polarization does not exeed 7%, while the optical degree of 
polarization >20% quite often; this is qualitatively consistent with the mm-wave polarization being 
affected by synchrotron self-absorption, which lowers the degree of polarization and rotates the 
EVPA by 90°; (ii) the magnetic field in the core region has a helical structure, with the optical and 
mm-wave emission regions located in different parts of the helix; and (iii) the jet has a spine-and- 
sheath structure such that, in a quiescent state, the optical polarization comes from a transition 
zone between the fast spine st ream and slower sheath flow (such emi ssion would be highly polarized 
perpendicular to the outflow; iLaing Ill980l : iD'Arcangelo et al.ll2009l ). while the mm-wave emission 
region is dominated by turbulence and has low polarization. In an active state, a disturbance in 
the jet defines the optical and mm-wave polarization properties: at millimeter wavelengths the 
disturbance orders the magnetic field perpendicular to the jet flow, while (/^opt depends on whether 
the disturbance or shear layer dominates the polarized flux. 



10. Discussion 

We have found a connection between mm-wave and optical outbursts and disturbances prop- 
agating down the radio jet in 3C 454.3. At least three major optical outbursts, Oki,Ok2, and 
Okh, observe d in 2 005-2008 are related to superluminal knots identified in the VLBA images. 



Jorstad et al.l (|2007l ) have found that in blazars the polarization properties of superluminal knots 
at 43 GHz are consistent with weak shocks and the optical polarized emission originates in shocks, 
most likely situated between the 86 GHz and 43 GHz VLBI cores. It appears that shocks and their 
interaction with underlying jet structure might be generally responsible for optical and high energy 
outbursts. 



10.1. Location of the Optical Emission in the Jet 

Comparison of the timing of optical outbursts and epochs of ejection of superluminal knots 
from the core can be used to determine the location of the sites where kinetic and internal en- 
ergy is dissipated in the form of radiation. The delay between the peak of an outburst and the 
time of the ejection of the associated superluminal component, along with the apparent speed 
of the knot, determines the location of the optical outburst in the jet with respect to the core, 
^'^opt = /5appcATopt/ sin0o. According to Tables [1] and HI the peaks of outbursts OkijOk2i and 
Ok3 occurred 12.6±5.6 pc upstream, 3.5±3.7 pc downstream, and 2±18 pc upstream of the core. 
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respectively. The average location for Ok 2 and Ok 3 is 2.8 ifc3.4 pc downtream of the core. In the 
case of a conical structure of the jet, iJorstad et al.l (|2005l ) have estimated the half opening angle 
of the jet in 3C 454.3 to be = 0.8° it 0.2°. We have obtained an average size of the core during 
our period of observations using modelling by circular Gaussians, acore=0.068±0.011 mas. This 
value of ttcore gives an upper limit for the transverse size of the core, which, along with the opening 
angle of the jet, provides an estimate of the distance of the mm- wave core from the central engine, 
Rbu <18±3 pc. Therefore, if the perturbation that created Kl is reponsible for Oki, then the 
maximum of the outburst took place when Kl was located at a distance 5.4±5.2 pc from the central 
engine, while the maxima of Ok2 and Ok3 outbursts occurred at a larger distance from the BH, 
21±4 pc, perhaps when knots K2 and K3 passed through the mm-wave core. This implies that 
the dissipation site for Oki was closer to the central engine (BH) than that for outbursts Ok2 and 
Oks- These results are consistent with the recently advanced ideas that the dissipation sites of 
optical and hi gh energy outbursts can occur at different l ocations and involve differen t sources of 
seed photons (jchisellini et allboOTl : iMarscher et alJbood . boiol : 



Chatter iee et al.l 120081 ). 



10.2. The Structure of the mm- Wave Core 

The most powerful optical outburst, Oki, preceded the passage of Kl through the 43 GHz 
core (ejection) by ~50 days. In addition, there are optical flares (RJD: 3537 and RJD: 3561, see 
Table[3]) that occurred within la uncertainty of the time of ejection of Kl. Outbursts Ok2 and Ok3 
were simultaneous (within la uncertainty) with the passage of K2 and K3, respectively, through 
the core. We isolate segments of the optical light curve within the time interval —aTo and +2aTo 
of each of knots Kl, K2, and K3 (Fig. [23l the values of aTo are given in Tabled]). Each ejection 
is accompanied by a series of flares (denoted in Fig. [23] as 1, 2, and 3) that starts within la of 
Tq. The duration of the series (time between the first and third peaks) is comparable, 50, 50, and 
65 days, for events Kl, K2, and K3, respectively. 

The profile of event Oki is very different from the profiles of Ok2 and Oks. The average of 
the rise and decay times for Oki has a width of ~47 days (Table [3]) , with a rise time of 68 days 
and a decay of only 26 days. Outbursts Ok2 and Oks, as well as the flares connected with the 
passage of Kl through the core, have much shorter durations and similar rise and decay times (see 
Fig. [23]) . A possible interpretation is that, in the case of Oki, located deeply within the ACZ, 
the energization of relativistic electrons is gradual and continuous, while for outbursts occurring 
in the VLBI core, relativistic electrons are energized abruptly owing to interaction between the 
disturbance and the core. The latter implies that the three-flare structure of optical outbursts seen 
in Figure [23] is related to the physical structure of the mm-wave core. Indeed, the super-resolved 
images of 3C 454.3 (Fig. [6]) support the idea that the core at 43 GHz might have a three-component 
structure. 



The core migh t be a syste r n of a lternat ing conical shocks and rarefactions, as suggested by 
Gomez et al. (jl997l ). Marscher ( 2006 ). and Cawthorne (j2006l ). In this case, the distance be- 
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tween constriction s , Zmax , can be calculated for an ultra-relativistic equation of state according to 
Daly Marscher I (jl988l 1: 

^max ~ 3.3 X r aJ^Qj.,,///, where fl*ore is transverse radius of the core in 



mas and rj = Pext/Pojet is the ratio of the externa l press ure to initial pressure (i.e., at the upstream 
boundary of the core) in the jet. iJorstad et al.l (|2005l ) have found that the relation between the 
viewing angle and Lorentz factor in blazars is consistent with r] ~3. If the core of 3C 454.3 is a 
system of three conical shocks, then the longitudinal size of the core in projection on the sky can be 
determined as acorcf""^ = 3 z^ax sin Bo- On the other hand, if the three-flare structure of the optical 
outbursts observed in the vicinity of the ejection is caused by a disturbance moving through this 



system of conical shocks, then a, 



',obs 



AT^fi, where ATg is the duration of the optical outbursts 



and n is the proper motion of the knot. Table [8] lists values of Ocorcf"'^ and aio^^ for events Kl, K2, 
and K3, which are consistent with the proposed association between a propagating disturbance 
and with the multi-component structure of the VLBI core providing a viable explanation for the 
three-flare pattern observed in the optical light curve. From standard gas dynamics, the presence 
of multiple standing shocks in three dimensions requires a high level of azimuthal symmetry of the 
pressure at the boundary of the jet in order for the structure to remain intact. 



10.3. Location of Millimeter- Wave Outbursts in the Jet 

We have derived the locations of 230 GHz outbursts Mki, Mk2, and Mks in the same manner 
as for optical outbursts Oki, Ok2, and Ok3 owing their association with superluminal knots. The 
maxima of Mki, -Mk2j and Mks occurred when a disturbance was located 4.6ib6.1 pc upstream, 
8.2ib4.8 pc downstream, and 15ib37 pc downstream of the core, respectively, with the average value 
3.8ib4.5 pc downstream of the 43 GHz core. This implies that the passage of a disturbance through 
the 230 GHz core should occur during the rising branch of a mm- wave outburst. 

The sharp mm-wave outbursts Ml and M2 (see Fig. [T]) seem unrelated to the appearance of 
new knots although we cannot reject the possibility that new knots were ejected near the time of Ml 
and M2 but could not be resolved on the images because they were weaker than knots Kl and K3, 
which were ejected before outbursts Ml and M2, respectively. However, the VLBA observations 
support the idea that the outbursts are caused by curvature in the inner jet that decreases the angle 
between a disturbance {Kl or K3) and the line of sight so th at the flux peaks wh en the disturbance 



passes through the minimum viewing angle, as suggested by lVillata et al.l ()2007l ). Indeed, Figure [9] 
shows curvature in the trajectories of Kl and K3 at a distance ~0. 05-0. 07 mas from the core. This 
is the location of knots Kl and K3 when outbursts Ml and M2 peaked. Comparison of Tables [T] 
and [3] indicates that Ml peaked at a time close to the maximum flux of knot Kl, and that the 
maximum of M2 occurred when knot K3 had the highest flux. The mm-wave spectral index during 
outbursts Ml and M2 corresponds to optically thin emission (see §8), as expected for radiation 
coming from downstream of the core. In addition, Table [3] shows that the timescale of Ml is ~4-5 
times shorter than li^mm of M2, consistent with Kl moving closer to the jet axis (a smaller radius 
of curvature) and K3 moving along periphery of the jet (a larger radius of curvature), which agrees 
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with the viewing angles of the knots (Table [T]). We use the delay between the times of ejection 
of Kl and peak of Ml, along with the kinematic parameters of Kl, to derive the location of the 
curvature of the jet, which is 54zbl6 pc downstream the core. 



10.4. Location of the High Energy Emission in the Jet 



The analysis of spectral characteristics of outbursts in 3C 454.3 reveals that neither during 
outbursts nor in quiescent states is the X-ray emission a continuation of the optical synchrotron 
spectrum. The most likely mechanism of X-ray production is inverse C o mpton scat tering of low 



energy photons by relativistic electrons, as suggested by ISikora et al.l (|2008l ) and iRaiteri et al 
(l2008bl ). 



Figure [2] displays normalized high-energy light curves, with the optical light curve superposed, 
during outburst Oki, the period covering Ok2 plus Oks, and M2. Table [2] indicates that the 
strongest correlation between optical and X-ray variations was during Oki, with a delay <1 day. 
Inspection of Figure [2] reveals that this correlation was dominated by the first ~15 days of the 
decaying branch of the outburst. (Unfortunately, the rising branch was not observed at X-rays 
and has a seasonal gap at optical frequencies; see Fig. [1]). According to the discussion above, the 
optical outburst took place upstream of the 43 GHz core. The X-ray outburst could have been 
pr oduced via the EC mechanism, with the seed photons coming from the dusty torus, as argued 
by ISikora et al.l (j2008l ) . Similarity in the decay at optical and X-ray frequencies implies that the 
primary cause of the outburst was an increase in the number of rel ativistic electrons. Ho wever, the 
EC model predicts that the X-ray spectrum should be rather flat (jGhisellini et al.l 120071 ). whereas 
ax is steepest during Oki (Table E]). Although this is comfortably less than aoptj it is striking that 
the X-ray spectrum was flatter during the quiescent state. This can be understood if the X-ray 
emission during the quiescent state is dominated by the EC process while the SSC mechanism 
contributes signiflcantly in the X-ray production during Oki • A possible delay of X-ray variations 
by_~l day relative to optical variations (Table [2]) supports such a hypotheses (jMcHardv et al 



19991 ). 



We have suggested that the optical flares observed within la uncertainty of the time of ejection 
ol Kl (Fig. [23|l originated in the VLBI core as the result of compression of the disturbance by 
standing shocks in the core. Figure [2] shows that two X-ray flares (designated as XI and X2) 
occurred during the same time span, which implies that these X-ray flares originated in the mm- 
wave core as well, probably via the SSC mechanism. The SSC model is supported by tentative 
delays of several days relative to the optical peaks for flares XI and X2, as can be inferred from 
Figure [2j The higher amplitude of the X-ray flares relative to the synchrotron flares (cf. Fig. [2]) is 
also expected in the SS C case if an increase in the number of relativistic electrons causes the events 
dChatteriee et ahlboosl ). 



The measurements during the Ok2+Ok3 period (middle panel of Fig. [2]) show an increase in 
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X-ray activity during Ok2 and possibly Ok3 ■ IVercellone et al.l (|2009l ) have found a good correlation 
between 7-ray and optical variations at the beginning of Ok3; and argue that the 7-ray emission 
is dominated by EC scattering of photons from the broad-line region by relativistic electrons in 
the jet. Although we agree that EC is probably the dominant mechanism for 7-ray production, 
we argue that the optical outburst Oks appears related to the 43 GHz core located ~20 pc from 
the broad-line region. We propose that seed photons for the 7-ray production are local to the 
mm-VLBI core, arising from a slower sheath surrounding the spine, although synchrotron radiation 
from th e disturbance itself pl us the mm-wave core should contribute to the production of 7-rays 
as well (jMarscher et al.ll20ld ). The presence of such a sheath i s supported by the wide range of 
apparent speeds observed in the jet, from 0.12 to 0.53 mas vr~^ (jJorstad et al.ll2005l ). 



The prominent mm-wave outburst M2 (see Fig. [T]) featured moderate 7-ray, X-ray, and optical 
flares. This high-frequency activity was especially pronounced during the early portion of the mm- 
wave outburst, but tapered off during its later stages. It is possible that a new superluminal knot 
was ejected during outburst M2, but could not be separated from K3 on the VLB A images. In 
this case, the multiple peaks of similar amplitude at 230 GHz, as well as the secondary optical 
maxima, can be explained by a knot passing through several standing shocks in the core region, as 
occurred during the passage of Kl, K2, and iC3 (Fig. [23]) . In the absence of a new ejection, the 
mm-wave outburst is probably not associated with the optical variability, rather it is caused by a 
change of the viewing angle of K'i when the knot was ~50 pc from the core (see §10.3). The optical 
emission could originate in the core according to polarization measurements, in a similar manner 
as proposed for the quasar 0420-014 (jD'Arcangelo et al.ll2007l ). 

The 7-ray light curve during event M2 (bottom panel of Fig. [2]) is well correlated with the 
variations in optical flux, and the X-ray variations correlate with both the optical and 7-ray light 
curves, albeit less strongly. The 7-ray and optical variations were simultaneous within a day 
(Table [2]), while the X-ray variations preceded the optical variations by 3ib2 days, although the 
X-ray measurements are sparse. The value of the X-ray spectral index during M2 is similar to ax 
during Oki- We suggest that the X-ray emission during M2 was produced via the SSC mechanism 
by the scattering of a wide spectrum of synchrotron photons — from far-IR to optical frequencies — 
by relativistic electrons with Lorentz factors 7~102"3. This can explain the moderate correlation 
between the X-ray and optical variations and smoother X-ray variability. If the delay of the optical 
with respect to the X-ray variations i s real, it can be understood in the manner proposed for the 
quasar 3C 279 (jChatterjee et al.ll2008l ). in which relativistic electrons are accelerated gradually. In 
this case, the acceleration mechanism should be different from that for Ok2 and Oks, and probably 
involves turbulence in the core. This is consistent with a significantly lower degree of polarization 
in the 43 GHz relative to optical polarization if the two emission regions are co-spatial but occupy 
different volumes (Table [6]). 

Our analysis of the 7 -ray and optica l light curves confirms a strong correlation between the 
two wavelengths found by iBonning et al.l (120091 ). The correlation persists b eyond M2, with 7-ray 
and optical variations having similar amplitudes of variability (see Fig. [2]). iBonning et al.l (|2009l ) 
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have concluded that the 7-ray outburst during M2 was dominated by EC scattering of IR/optical 
photons from the accretion disk and/or broad-hne region by electrons with 7 ^ lO^-"^. Although this 
may seem a reasonable explanation for the observed strong correlation between 7-ray and optical 
variations, polarization observations (see Fig. [3] and Table [6]) indicate that during M2 (RJD: 4550- 
4840) the position angle of the optical polarization tended to align with (^943 in the core. This 
implies that the variable optical emission arose in the vicinity of the core, as discussed above. For 
this reason, we favor the alternative scenario, wherein the 7-rays are produced by scattering of 
synchrotron seed photons from the sheath of the jet. 

We can be more specific regarding the natu re of the seed photo ns by considering the steepening 
of the 7-ray spectrum above 2 GeV reported bv lAbdo et al.l (|2009bl ) . from a spectral index of 1.3 to 
2.5, along with steepening of the spectral index of synchrotron emission from 1.4 at near-IR to 2.3 
at soft UV wavelengths (Table [5]). We can explain the steepening of the 7-ray spectrum, as well as 
the rapid fall-off in the synchrotron emission toward higher ultraviolet frequencies by a steepening 
of the electron energy distribution above an energy 7 ~ lO'* in rest- mass units , where we assume 
a magnetic field strength < 1 G. Recent theoretical calculations by Reynolds (j2009l ) show that, 
for an inhomogeneous source, more pronounced steepening than by 0.5 in the source's integrated 
spectral index is possible due to a combination of synchrotron losses and geometrical effects. If the 
2 GeV 7-rays are produced by the EC process, the seed photons should hav e observed frequen cies 
near 10i2(r/25)-2((5/25)-2 Hz. This favors either the putative dust torus dSikora et ahlboosl ^ or 
sheath of the jet as the source of the seed photons. The similarity of the 7-ray and optical light 
curves, including very rapid variations (see §7), can then be explained as a consequence of the 
variations being caused by changes in the number of electrons with 7 ~ 3 X lO^-lO^ 

over a limited 

volume within a given disturbance. 



11. Conclusions 



The optical. X-ray, and radio lig ht curves of blazars that have beeri assembled and ana- 
lyzed over the past several years (e . g.. IVillata et al.ll2004l : iBach et al.l l2006l : iRaiteri et al.l l2008bl : 



Chatteriee et al.1 12OO8I : IVillata et all l2009l : iMarscher et al 



20101 ) strongly support the conclusion 



that the long-term variability at all wavelengths is governed by physical activity in the relativistic 
jet that we can image with the VLBA. Our multifrequency study of the quasar 3C 454.3 during 
the dramatic activity of 2005-2008, including intense monitoring of the structure of the parsec-scale 
jet, reveals strong connections among the observed optical, X-ray, and 7-ray outbursts, its spectral 
and polarization behavior, and disturbances propagating down the jet. 

We have found that the location of the most prominent X-ray and optical outburst (in 2005) 
was within~10 pc from the BH. The event was the result of a disturbance propagating down the jet 
through a photon field that included synchrotror i radiation from b oth the high-P spine and slower 
sheath of the jet, and perhaps a hot dust torus (ISikora et al.ll2008l ) with X-ray emission produced 
by both the SSC and EC mechanisms. 
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Our intensive monitoring with the VLBA at 43 GHz reveals that interaction between distur- 
bances in the jet and the mm-wave core results in X-ray and optical outbursts and can cause optical 
intra-day variab ility. This implies that the mm-wave core is a physical feature of the jet, e.g., a 
standing shock (jCawthorne I l2006l : iJorstad et al.l 120071 : iD'Arcangelo et al.l 120071 ) , that is different 
from the cm-wave core, the location of which is probably determined by the jet opacity. 

We have observed three events when the passage of a superluminal knot through the core 
coincided with a series of optical outbursts and an increase in the X-ray emission. We infer that the 
structure of the optical outbursts reflects the fine-scale structure in the core, a scenario supported 
by super-resolved VLBA images. The X-ray production in the core region can be dominated by 
either the EC or SSC mechanism , depending on the properties of the disturbance, as suggested by 
Katarzynski fc Ghiselliiin koO'^ ). 



Vercellone et al.l (|2009l ) have reported a good correlation between 7-ray and optical variations 
in autumn 2007, which coincides with a series of optical outbursts that we relate to an event in the 
core. This suggests that 7-ray emission is also produced during interactions within the core. If the 
7-rays originate via the EC mechanism, then the seed photons should be local to the 43 GHz core 
and could come from the sheath of the jet. 

We have observed three episodes of rotation of the position angle of optical polarization near 
the times of passage of superluminal components through the core. We interpret this as evidence 
that the mm-wave core of 3C 454.3 is located near the end of the acceleration and collimation zone 



of the jet flow, where the magnetic field has a toroidal structure (e.g.. lVlahakis 1120061). Our findings 



agree with other studies reporting such rotations related to disturb ances in the jet ([Marscher et al 



20081 : iLarionov et al.ll2008l : iMarscher et al.ll20ld : lAbdo et al.ll20ld ). This implies that rotation of 
the optical polarization position angle is a common occurrence during the early stages of flares in 
blazars. In the case of 3C 454.3, we tentatively infer an inverse relationship between the Lorentz 
factor of the superluminal component and the rate of rotation of the optical EVPA, as expected 
if higher values of F are reached farther from the BH where the jet has expanded to a larger 
cross-sectional radius. 

Our study of the 7-ray bright quasar 3C 454.3 demonstrates an obvious connection between 
events at different wavelengths, the complexity of their relationships, and the primary role of the 
parsec-scale jet in the generation of these observed events. The ability of the Fermi Large Area 
Telescope to produce well-sampled 7-ray light curves for the first time, combined with the ultra-high 
resolution imaging of the VLBA and flux and polarization curves at other wave bands, provides 
an unprecedented opportunity to locate the events within the most compact regions of the jets of 
blazars. This endeavor makes possible a deeper understanding of the dynamics of relativistic jets 
and physics of the radiative processes. 
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Tabic 1. Parameters of Moving Components 



Knot To To Smax Ts^^^ n fiapp r 5 Go 

yr RJD Jy RJD mas yr"^ 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



Kl 2005.50±0.08 3553±29 7.00±0.34 3772 0.10±0.02 4.9±0.9 12.3 23.4 1.0 
K2 2007.49±0.07 4279±25 3.81±0.56 4343 0.18±0.05 8.3±2.3 10.7 17.6 2.5 
K3 2007.93±0.10 4439±35 11.74±0.52 4720 0.09±0.03 4.1±1.4 24.7 49.1 0.2 



Note. — Columns: 1 - name of component; 2 - time of ejection from the core in years; 3 - time of ejection from 
the core in RJD; 4 - maximum fiux of component, 5 - epoch of the maximum flux of component; 6 - proper motion; 
7 - apparent speed; 8 - Lorentz factor; 9 - Doppler factor; 10 - viewing angle 



Table 2. Results of Correlation Analysis 



RJD:3502-3619 RJD:4292-4450 RJD:4613-4858 



Waves 


fmax 


T 


y^max 


T 




T 






days 




days 




days 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


7-ray/X-ray 










+0.75 


+2±3 


7-ray/Opt 










+0.91 


+0±1 


X-ray/Opt 


+0.92 


+1+1 






+0.48 


-3+2 


X-ray/230 


+0.88 


-42+3 






+0.58 


-16+3 


Opt/230 


+0.98 


-48±5 


+0.46 


-21±2 


+0.59 


-13±4 



Note. — Columns: 1 - two bands used for correlation analysis; /max - maximum coefficient of linear correlation; 
T - time delay between the two bands at the maximum coefficient, negative delay corresponds to higher frequency 
variations leading 
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Table 3. Parameters of Total Flux Outbursts 







Optical 






230 GHz 






37 GHz 




Name 


Topt 


Eopt 


Wopt 


2230 


-E230 


W230 


T37 


Esr 


W37 




RJD 


mjy day 


days 


RJD 


Jy day 


days 


RJD 


Jy day 


days 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Oki/Mki 


3501±1 


2194 


47 


3534±5 


9343 


104 


3700±5 


6200 


180 




3537±1 


166 


15 
















3561±1 


115 


18 
















3673±1 


96 


19 














Ml 








3781±5 


1140 


38 


3788±5 


528 


24 


Ok2/Mk2 


4301±1 


373 


19 


4330±5 


1985 


50 


4331±5 


1871 


104 




4334±1 


170 


12 
















4382±1 


113 


16 














Ok3/Mk3 


4437±1 


410 


16 


4454±5 


1020 


51 


4450±5 


2090 


64 




4491±1 


253 


27 
















4661±1 


164 


16 














M2 








4705±5 


6810 


134 


4707±5 


6663 


155 



Note. — Columns: 1 - designation of outburst; 2, 3, 4 - epoch of the peak, power, and time scale of optical 
outburst, respectively; 5, 6, 7 - epoch of the peak, power, and time scale of outburst at 230 GHz, respectively; 8, 
9, 10 - epoch of the peak, power, and time scale of outburst at 37 GHz, respectively; uncertainties of the epochs of 
peaks correspond to the smoothing time of light curves 



Table 4. Connection between Ejection and Peak of Outbursts 



Knot 
(1) 



ATopt 
days 
(2) 



AT230 

days 
(3) 



Ar37 

days 
(4) 



Kl -52±30 -19±34 147±34 
K2 22±26 51±30 52±30 

K3 -2±36 15±40 11±40 



Note. — Columns: 1 - name of component; 2, 3, 4 - delay of the peak of optical, 230 GHz, and 37 GHz outbust, 
respectively, with respect to the ejection time of component, AT^, =Ti, — To 
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Table 5. Optical, near-IR, UV, and X-ray Spectral Indices" of the Synchrotron 

Emission Components 



State 
(1) 


(2) 


"opt 

(3) 


syii 

(4) 


"UV2 

(5) 


(6) 


Epochx (RJD) 
(7) 




1.88±0.03 


1.67±0.11 






0.72±0.04 


3505.5814 


Ok2+Ox3 


1.31±0.02 


1.68±0.01 


2.31±0.14 


0.81±0.17 






Ok2 










0.41±0.13 


4299.1418 


Ok3 










0.56±0.12 


4426.7780 


M2 


1.65±0.02 


1.64±0.01 


2.33±0.22 


1.16±0.15 


0.66±0.18 


4706.2752 


QS" 


1.42±0.01 


1.97±0.03 






0.57±0.01 


4087.5 



Note. — Columns: 1 - state of activity; 2 - IR spectral index of a synchrotron component; 3 - optical spectral index 
of a synchrotron component; 4 - UV (2500-3500 A) spectral index of a synchrotron component; 5 - UV (2000-2500 
A) spectral index of a synchrotron component; 6 - X-ray (2.4-10 keV) spectral index; 7 - epoch of measurement 

''AH spectral indices are defined as Sv oc v^" 

The values of Oopt and Oir are from lHagen-Thorn et al. I l|2009l ) 



'^The values of are from lRaiteri et al.l (|2007l l 
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Table 6. Quasi-Simultaneous Optical and 43 GHz Polarization Measurements 





Popt 


V^opt 




p43 


b 

¥^43 


Pi- t 


RJD 


/O 


deg 


RJD 


/o 


deg 


deg 




(2) 




(4) 


C5~l 




(7) 


*3563.4871 


4.90±0.20 


-134.4±1.1 


3568 


1.0±0.2 


-127.0±4.4 


-94.8±3.4 


3656.3266 


10.72±1.45 


-145.9±3.9 


3650 


1.7±0.3 


-95±4.9 


-98.4±4.4 


3684.2745 


8.09±2.06 


-21.3±7.3 


3689 


0.8±0.2 


7.1±4.6 


-89.0±5.3 


3697.3281 


9.89±1.99 


-26.1±5.8 


3695 


2.1±0.4 


9.5±2.9 


-101.2±7.5 


3729.2215 


8.93±0.58 


17.1±1.9 


3726 


2.2±0.2 


-0.5±2.4 


-86.6±3.3 


* 3906.4829 


6.46±3.50 


-114.6±15.6 


3908 


1.8±0.3 


-136.7±5.5 


-90.9±2.9 


3932.5166 


12.10±6.53 


-59.3±15.5 


3935 


1.6±0.5 


-92.3±8.4 


-90.6±2.0 


*4014.3535 


15.83±6.40 


-30.5±11.6 


4014 


1.1±0.5 


-22.5±14.4 


-101.0±4.6 


* 4050. 7222 


2.00±0.70 


-53.3±10.0 


4057 


1.4±0.5 


-49.0±10.7 


-100.4±5.3 


4086.1988 


9.31±0.82 


-161.9±2.5 


4087 


0.6±0.5 


-67.9±22.8 


-94.7±2.1 


*4269.4898 


2.66±1.09 


-124.1±11.8 


4266 


0.7±0.2 


-126.0±5.9 


-86.7±6.5 


*4300.4824 


21.44±0.24 


-100.8±0.6 


4294 


2.6±0.3 


-97.0±3.4 


-86.9±4.7 


4319.4380 


5.76±0.61 


-58.2±3.0 


4320 


3.2±0.4 


-127.0±3.4 


-103.6±3.3 


4343.4985 


7.05±1.85 


-61.7±7.5 


4343 


3.1±0.3 


-111.8±3.0 


-102.2±4.1 


*4405.2629 


6.48±0.53 


-117.0±2.4 


4406 


1.2±0.3 


-120.7±7.4 


-103.6±2.8 


4627.4868 


7.82±0.55 


-12.9±2.0 


4629 


2.0±0.5 


-129.7±6.5 


-90.3±3.4 


4654.4839 


8.61±0.52 


-147.8±1.7 


4653 


1.4±0.4 


-109.0±7.4 


-104.4±2.2 


4794.6152 


2.60±0.07 


-143.1±0.7 


4787 


1.4±0.5 


-155.3±9.5 


-123.2±3.0 


*4828.6533 


0.79±0.16 


-155.4±4.1 


4822 


2.5±0.3 


-162.0±3.9 


-112.6±2.5 



Note. — Columns: 1 - epoch of optical observation; 2 - optical degree of polarization; 3 - optical position angle 
of polarization; 4 - epoch of VLBA observation; 5 - degree of polarization in the core region at 43 GHz; 6 - position 
angle of polarization in the core region at 43 GHz; 7 - direction of the inner jet in projection on the sky plane 

''symbol (*) marks epoch at which the optical polarization position angle and EVPA in the 43 GHz core align 
within 1 a uncertainty 

''position angles of polarization at 43 GHz are corrected for RM according to lJorstad et al.l ( 20071 1 
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Table 7. Quasi-Simultaneous Optical and 86 GHz Polarization Measurements 



Epochopt 


Popt 


V^opt 


Epochse 


P86 




ejet 


RJD 


% 


deg. 


RJD 


% 


deg. 


deg. 


Uj 






(A\ 






\l) 


4304.5098 


13.84±0.94 


-82.1±1.8 


4304 


2.61±0.59 


-110.4±6.5 


-86.9±4.7 


4313.4829 


6.18±0.83 


-78.0±3.9 


4314 


6.93±1.19 


-112.2±4.9 


-103.6±3.3 


*4342.4062 


3.70±1.31 


-110.9±10.1 


4339 


4.09±0.54 


-103.4±3.8 


-102.2±4.1 


4346.4395 


9.35±1.37 


-28.4±4.2 


4346 


3.62±0.57 


-105.9±4.5 


-102.2±4.1 


4355.4492 


3.57±1.93 


-3.5±9.7 


4353 


2.80±0.47 


-96.8±4.8 


-102.2±4.1 


4359.4316 


12.91±3.87 


-41.9±6.6 


4357 


2.11±0.57 


-100.0±7.7 


-103.5±2.0 


4371.3677 


26.98±1.89 


-51.6±2.0 


4372 


1.86±0.56 


-91.5±8.6 


-103.5±2.0 


4395.3691 


23.48±4.11 


-169.1±5.0 


4393 


1.45±0.46 


-102.4±9.2 


-103.6±2.8 


*4410.3267 


14.09±1.01 


-99.5±2.1 


4408 


1.46±0.47 


-87.3±9.2 


-103.6±2.8 


4768.4941 


8.26±1.72 


-61.0±8.6 


4768 


1.39±1.02 


-7.9±21.1 


-132.2±3.5 


4769.5107 


4.11±0.60 


-58.0±6.0 


4769 


1.79±0.47 


-169.6±7.5 


-132.2±3.5 


4773.7139 


5.20±1.42 


-79.4±7.8 


4773 


1.46±0.47 


-167.1±9.1 


-132.2±3.5 



Note. — Columns: 1 - epoch of optical observation; 2 - optical degree of polarization; 3 - optical position angle 
of polarization; 4 - epoch of observation at 86 GHz; 5 - degree of polarization at 86 GHz; 6 - position angle of 
polarization at 86 GHz; 7 - direction of the inner jet in projection on the sky plane 

^symbol (*) marks epoch at which the optical polarization position angle and EVPA at 86 GHz align within 1 a 
uncertainty 



Table 8. Longitudinal Size of the VLBI Core 



Knot 

(1) 


2!max 

mas 
(2) 


Ijtheor 
"•core 

mas 
(3) 


AT, 

days 
(4) 


Ijobs 
t*core 

mas 
(5) 


Kl 


0.50 


0.026 


50 


0.013 


K2 


0.56 


0.041 


50 


0.025 


K3 


1.01 


0.011 


65 


0.018 



Note. — Columns: 1 - name of event; 2 - distance between constrictions; 3 - estimate of longitudinal size of the 
core from the theoretical model; 4 - estimate of longitudinal size of the core from observational consideration 
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Fig. 1. — Light curves of the quasar 3C 454.3 at different frequencies. Dotted hues show times 
of ejections (times of coincidence with the position of the VLBI core at 43 GHz) of superluminal 

components Kl, K2, and K3 and sohd hne segments show uncertainties in the ejection times; 
second panel: open circles - Liverpool telescope data in [/-band, open triangles -Swift UVOT data 
in [/-band, and filled triangles - Swift UVOT data in UWl-hand; 345 GHz (open circles, fourth 
panel) and 43 GHz measurements (open circles, fifth panel) are shifted by +15 Jy for clearity. 
Symbols designate the most prominent outbursts at different wavelengths. 
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Fig. 2. — Normalized i?-band (black), X-ray (green, i?XT£' measurement - circles and Swift mea- 
surements - triangles), and 7-ray (blue) light curves for three periods. The maximum flux within 
each panel is adopted as the normalization factor. Top panel: S'J^'^^=63.6 mJy, £'^^'^^=10. 13 /uJy 
(X-ray flux is at 4 keV); middle panel: 5^^^=37.17 mJy, 5^^^=2.138 cts/s (X-ray flux is at 0.3- 
10 keV); bottom panel: 5^'''^=17.82 mJy, 5^^^=1.614 cts/s, 5;^''^=1.2x 10"^ phot/cmVs (X-ray 
flux is at 0.3-10 keV , 7-ray flux i s at 0. 1-300 GeV). Blue triangles show the 7-ray measurements 
taken from Fig. 1 of lAbdo et al.l ()2009bl ). Error bars for the optical observations are less than the 
size of symbols. 



- 38 - 



50 



1 □□ 



^ 2D 

^ in 

□ 

3D 



2D 




- n— band 



ID - 



8 
6 



IDD 
D 

-lm 



- Jtn 



D 1 1 h 



: 43 QHz Cqi^BE QHz 



H 1 1 h 



ilk 



-t—=] — ^ 



H 1 1 h 



7 R-tmrri 

= 43QHzGBrq^BE QHz 



H 1 1 1- 




H — I — H — I- 



H 1 h 





X3DD 



35DD 4DDD 
JD 245]D[X>4- 



Fig. 3. — Flux and polarization curves of the quasar 3C 454.3 at different frequencies: top panel 
- total flux density in R-band; second panel - degree of polarization at optical wavelengths; third 
panel - degree of polarization at 86 GHz (green circles) and in the VLBI core at 43 GHz (black 
circles); and bottom panel - position angle of polarization at optical wavelengths (black circles), at 
86 GHz (green circles), and in 43 GHz VLBI core (red circles); dotted lines show times of ejections 
of superluminal components and solid line se gments show uncerta inties in the ejection times. Values 
of ¥57111111 are corrected for RM according to IJorstad et al.l ()2007l ) . 
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Fig. 4. — 43 GHz total intensity images of the quasar 3C 454.3 with S'pcak=19.43 Jy/beam and 
beam=0. 14x0.30 mas^ at PA=-W°. Contours represent 0.0625, 0.125, 0.25,. ..64% intensity of the 
peak. 
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Fig. 5. — 43 GHz total (contours) and linearly polarized (color scale) intensity images of 3C 454.3 
with different convolving beams; top : 0.33x0.14 mas^ at PA=-10°, S'peak=6.20 Jy, 5'Pg^k=0.076 Jy; 
bottom: 0.10x0.10 mas^, S'peak=5-66 Jy, S^^^]^=Q.QQ'l Jy. Line segments show the plane of polar- 
ization and letters identify components in the jet. 
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Fig. 6. — 43 GHz total (contours) intensity images of 3C 454.3 with convolving beam 
0.05x0.05 mas^ at different epochs; top : 2009 January 24 (RJD:4856) 5peak=5.10 Jy, 

S'P^^j^=0.067 Jy; bottom: 2007 November 1 (RJD:4406) S'peak=5.84 Jy 5'P^,^j^=0.050 Jy Line seg- 
ments show the plane of polarization and letters identify components in the jet. The tripod of lines 
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Fig. 7. — Distance of components from the 43 GHz core as a function of time. The size of each 
symbol is proportional to the logarithm of the flux density of the knot. Dotted lines indicate an 
approximation for the presumably ballistic motion of components K l, K2, and K3 near the core. 
The dashed line shows the position of stationary knot C reported by lJorstad et al.l 1)20051 ) . 
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Fig. 8. — Direction of the jet within 0.3 mas from the core as a function of time. Dotted hnes show 
the times of ejections of superluminal components Kl, K2, and K3 from the 43 GHz core and solid 
line segments show uncertainties in the ejection times. 
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Fig. 9. — Light curves (top) and trajectories {bottom) of moving knots; the cross indicates typical 
uncertainties in position of components {upper left coner of the bottom panel) 
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Fig. 10. — 43 GHz total (left) and polarized (right) intensity images of the quasar 3C 454.3 
during period when Kl emerged from the core, 5pcak=8-25 Jy/beam, 'S'peak"^-^^-'- Jy/beam, 
beam=0. 1x0.1 mas^. Total intensity contours at 0.25,0. 5, ...64 % of the peak and polarized in- 
tensity contours at 2,4,... 64 % of the peak. Sticks over the polarized intensity contours indicate the 
plane of polarization. 
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Fig. 11. 43 GHz total (left) and polarized (right) intensity images of the quasar 3C 454.3 dur- 
ing period when K2 emerged from the core, S'pcak=15.69 Jy/beam, S^^^^=0.2Q8 Jy/beam. Total 
intensity contours at 0.25,0. 5, ...64 % of the peak and polarized intensity contours at 2,4, ...64 % of 
the peak. Sticks over the polarized intensity contours indicate the plane of polarization. 
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Fig. 12. — 43 GHz total (left) and polarized (right) intensity images of the quasar 3C 454.3 
during period of K3 emerging from the core, 5peak=14.28 Jy/beam, "Speak "^•■'-^^ Jy/beam, 
beam=0. 1x0.1 mas^. Total intensity contours at 0.25,0. 5, ...64 % of the peak and polarized in- 
tensity contours at 2,4,... 64 % of the peak. Sticks over the polarized intensity contours indicate the 
plane of polarization. 
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Fig. 13. — Cross-correlation function between X-ray, optical, and 230 GHz light curves for the 
indicated time intervals. A negative delay corresponds to higher frequency variations leading. 
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Fig. 14. — Left: Cross-correlation function between X-ray, optical, and 230 GHz light curves for 
the time interval RJD: 4610-4860. Right: Cross-correlation function between 7/X-ray variations for 
the time interval RJD: 4689-4870, and between 7-ray and optical light curves for the time interval 
RJD: 4650-4870. A negative delay corresponds to higher frequency variations leading. 
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Fig. 15. — Cross-correlation function between the 230 GHz and 37 GHz light curves for periods 
RJD:3000-4850 (solid line) and RJD:4000-4850 (dotted line). 
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Fig. 16. — Values of standard deviation of R-band magnitudes averaged within a day (filled circles). 
The average uncertainties of individual measurements within a day are shown by crosses and are 
given negative values to distinguish them from the standard deviations. Dotted lines show times 
of ejections of superluminal components. 
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Fig. 17. — Spectral index at 230-37 GHz (1-8 mm). Dashed line shows amm during a quiescent 
state, QS. Dotted lines show times of ejections of superluminal components and solid line segments 
show uncertainties in the ejection times. 




Fig. 18.— Flux-Flux dependences for UVOT bands for periods RJD: 4292-4450 (left) and JD 
4613-4858 (right). 
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Fig. 19. — Near-IR (left), optical {middle), and UV (right) relative spectral energy distribu- 
tions of synchrotron emission from 3C 454.3, normalized to K-band (Lg(z^o)=14.140), R-band 
(Lg(z^o)=14.670), and UWl-band (Lg(z^o)=15.056), respectively, during different events. 
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Fig. 20.— 



Spectral index at 0.1-300 GeV. 
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Fig. 21. — Distributions of position angle of polarization with respect to the jet direction at different 
wavelengths for all {top panels) and high polarizat i on (bo ttom panels) measurements; values of ipjmm 
are corrected for RM according to IJorstad et al.l ()2007l ) . 
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Fig. 22. — Polarimetric behavior at optical wavelengths (filled circles) and in 43 GHz core (open 
triangles) near the time of ejection of superluminal components (dotted lines), solid line segments 
show uncertainties in the ejection times. 
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Fig. 23. — Optical photometric variability within —laTo to +2(7 To time interval relative to the time 
of the ejection, To of each superluminal component (dotted line); values of To and their uncertainties 
are given in Table [H 



